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In order to discuss the well-posed initial value formulation of the teleparallel gravity and apply 
it to numerical relativity a symmetric hyperbolic system in the self-dual teleparallel gravity which 
is equivalent to the Ashtekar formulation is posed. This system is different from the ones in other 
works by that the reality condition of the spatial metric is included in the symmetric hyperbolicity 
and then is no longer an independent condition. In addition the constraint equations of this system 
are rather simpler than the ones in other works. 
PACS numbers: 04.20.Cv, 04.20. Ex, 04. 25.Dm 

I. INTRODUCTION 

As the closest alternative to general relativity (GR), teleparallel gravity can be traced back to Einstein [1] who 
regarded it as a unified field theory and attempted to use it to supersede GR. Teleparallel gravity can be regarded as a 
translational gauge theory [2-4] , which makes it possible to unify gravity with other kinds of interactions in the gauge 
theory framework. Poincare gauge theory is a natural extension of the gauge principle to spacetime symmetry , and 
represents a alternative to GR (for more general attempts see [3]). In particular, teleparallel gravity was regarded as a 
promising alternative to GR until the work of Kopczynski [5] , who found a hidden gauge symmetry which prevents the 
torsion from being completely determined by the field equations, and concluded that the theory is inconsistent. Nester 
[6] improved the arguments by showing that the unpredictable behavior of torsion occurs only for some very special 
solutions (see also [7]). Hecht et al. [8] investigated the initial value problem of teleparallel gravity and conclude that 
it is well defined if the undetermined velocity are dropped out from the set of dynamic velocities. 

Teleparallel gravity possesses many salient features. Because of its simplicity and transparency teleparallel gravity 
seems to be much more appropriate than GR to deal with the problem of gravitational energy momentum [4, 9, 10]. 
Nester [9] succeeded in proving the positivity of total energy for Einstein's theory in terms of teleparallel geometry. 
He found that special gauge features of teleparallel gravity, which are usually considered to be problematic, are quite 
beneficial for this purpose. Mielke [10] used the teleparallel approach to give a transparent description of Ashtekar's 
new variables [11]. Andrade et al. [12] formulated a five-dimensional equivalent of Kaluza-Klein theory. Although 
quantum properties of the Poincare gauge theory are, in general, not so attractive, the related behavior in the specific 
case of the teleparallel theory might be better [13], and should be further explored. 

The canonical Hamiltonian approach is the best way to clarify both the nature of somewhat mysterious extra gauge 
symmetries and the question of consistency of teleparallel gravity. It is found [7] that the presence of nondynamical 
torsion components is not a sign of an inconsistency, but a consequence of the constraint structure of the theory. 
In some works[7,14, 15], the Hamiltonian formulation of teleparallel gravity has been developed. However, the well- 
posed initial value formulation of teleparallel gravity has not been discussed as yet. As is well known, hyperbolic 
formulation of the Einstein equation is one of the main research areas in GR [16]. This formulation is used in the 
proof of the existence, uniqueness, and stability of the solutions of the Einstein equation by analytical methods. Thus 
far, several first-order hyperbolic formulations are proposed [17-20]. The recent interest in hyperbolic formulation 
arises from its application to numerical relativity [21]. It is proved that the Einstein equation in Ashtekar's variables 
constitutes a symmetric hyperbolic system [22, 23]. A question naturally arises whether there is a well-posed initial 
value formulation of teleparallel gravity which is equivalent to or different from the hyperbolic formulation of the 
Einstein equation. If it exists, can it give us some new perspectives and be applied to numerical relativity? An 
answer to this question will be given in this paper. A self-dual teleparallel formulation of general relativity which is 
equivalent to Ashtekar's formulation has been developed, its canonical Hamiltonian analysis has been given and used 
to clarify the gauge structure of the theory [15]. In this paper a new symmetrical hyperbolic system of the Einstein 
equation will be posed in terms of the two-spinor formulation based on [15]. A new fact we find is that the reality 
condition of the spatial metric is included in the symmetric hyperbolicity and then is not independent. In addition, the 
constraint conditions take a rather simple form. All of these reduces the number of independent conditions imposed 
on the equations and then simplifies the relevant problems largely. In Sec. II, a canonical formulation of the self-dual 
teleparallel equivalent of GR is given, the Hamiltonian (evolution) equations and the constraint equations are written 
in terms of the two-spinor dyad and the canonical conjugate momenta. In Sec. Ill, by introducing a new variable the 
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evolution equations are rewritten as first-order forms. Then in Sec. IV, the conditions of the symmetric hyperbolicity 
of the evolution equations are discussed and the relations between the conditions of the symmetric hyperbolicity and 
the reality conditions of the spatial metric are obtained. Section V is devoted to some conclusions.. 



II. THE CANONICAL FORMULATION OF THE SELF-DUAL TELEPARALLEL EQUIVALENT OF GR 

We start with the Lagragian of the self-dual (or chiral) teleparallel formulation of general relativity [15,10] which 
is equivalent to the Ashitekar Lagrangian [11] and is written in terms of the two-spinor formulation [24] 

V+ = Na[u)(AB) AC u DB dc - u {AB)cd uj cbad - V2u ± cdu {CE) e D ], (1) 

where ujab cd is the self-dual spin connection, uj(ab) CD is the Ashtekar variable, 

Uj_CD = n AB LOABCD = n AB LO[AB]CD, (2) 

and 

n AB = j_ {t AB_ N A B)) (3) 

with the determinant of the inverse soldering form a = det o-^ AB = jj\f—g , the lapse TV and the shift N AB . 
In the two-spinor formalism [15] the basic variable is chosen to be the dyad CaA- 

CoA = OA , ClA = M , 

and the self-dual spin connection luabcd can be expressed by 

WABCD = Cc b dAB(bD- (4) 

Using (2) (3) and (4) one gets 

W -LOD = ~J^^C b CbD -Jj NAB U(AB)CD, (5) 

where 

(bD = t AB 9AB(bD- 

Then the Lagrangian becomes 

V+ = Na[^A B) AC ^ D c ^ ab) cd^ CB)AD ] 

-2V2a( c b CbD ^ EC) e D + 2V2aN AB cj (A B)CD^ EC) E D . 
The canonical momentum conjugate to CbB is then 

Qy+ 

f D = — = -2V2aCc"^ EC) E D , (6) 
dCbD 

which leads to 

cb B 
W {AC) — -^UaV ■ 

The gravitational Hamiltonian can be computed 

H+ = $> D (to -V+= NH±_ + N ab Hab, (7) 

where 
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H± = ^e ab F C tfc + vCc<; bA d {AB) ( aD d( CB \ b D , (8) 

and 

Hab = d {AB) C aD p aD ■ (9) 
Computing the variation, we obtain the Hamiltonian equations 

<»A = p[ = + N CB d {CB) ( aA , (10) 



P 



KaA 

N CB d (BC) p aA + 2NaCDC bB d (BC) d {DC) C b A - ^C A P* B p bB 



^NaC A CcC bE d (EB) C cD d^C b D + SNaC D C bB C E d {BC) C C Ed^C b A 
~2NaC bc d {CB) C D d^ AB \ bD + 2Na( bB d {BC) C D d( DC \ b A 

+2NaCDd {BC) C hB d {DC) C b A + d (BC) N CB p aA - 2d (BC) NaC D C bB d^C b A 1 (11) 



and the constraint equations 

H±=0,Hab = 0. (12) 
The detail of the constraint structure of the theory can be found in [15]. 



III. THE FIRST-ORDER EVOLUTION EQUATIONS 

Since the Hamiltonian equations of p aA include the second-order terms 2Na( a D( bB d( B c)d( DC ) ( b A , in order to get 
the first-order evolution equations we decompose U(ac) DB hito its trace-free and trace part 

DB DB , c ~aB /iq\ 

W (AC) = ^trf(AC) + 4 ^/2 cr £c ^ aAP ( L6 > 

and introduce a new variable q AB and a real constant spatial 1-form ipcD by 

\/2 

Utrf(AC)DB = -g^1pD(CQA)B, (14) 

with the properties 

^CD = 1p(CD), QAB = Q[AB] • (15) 

Substituting (14) into (13) and using (4) we obtain 

d{AC)CaB = ^bPa(CQA)B ~ (o,(CPA)b\- (16) 



Then we can compute 



and 



■ N„ V2N^ _ V2N^ B _ 

(aA= -^PaA g~ QaBPCA H ^ VaBQCA, (17) 



a=-[iVC 0>1 p„A + ; ^— ^b«ca]. (18) 



Using these results we can rewrite (11) as 
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~aA _ \/2iV B c ^ A A jCBc, ~aA 



P = g C Cb 0(BC)P + M d{BC)P 

C w c*(bc)3 H 5— V 0(bc)Q 
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+Q P (p,q), (19) 

where source term Q p (p, q) is a quadratic polynomial of p aj4 and . 
Introducing a constant spatial vector <p AB satisfying 

^ACf AB = e c B (20) 

and using (16) one can obtain 

qAB = ^<j( bD ^D C d {AC) C b B - \vAcCa C p a B- (21) 

Taking the time derivative leads to 

8\/2 j-hj~) n 1 1 _> c — 

QAB = —^~ a C (AC ) CbB -^fAcCa P B 

, 8^2 _ bD c 8\/2 ■ bD c 

H 17 C <J(AC)QbB H — CT C, </5D 0(AC)ibB 

— VAC (a C p a B- (22) 

Using (17), (18), and (21) we get the evolution equation of qAB- 



Qab 



2 -^C D VD C d {AC) p aB ^c D VA c d iCD) p aB + \N DE VE c c D d {AC) p aB + l -N CE VAD c D d (CE) p aB 



and 



+Q,(p,g), (23) 

where source term Q g (p, q) is another quadratic polynomial of p aA and q AB . The source terms Q p (p, q) and g) 
do not contain any derivatives of the fundamental variables other than the lapse N and the shift N AB . 
Substituting (16) into (8) and (9) one can obtain: 

v% 

Wab = ~rz [ipcAqBD + iPcbciad]? 00 , 

l0(7 

H ± = -^PABP AB - ^J A BPACq BC ~ - 

And then the constraint equations (12) leads to 

qAB = -qBA, (24) 

and 

39p AB p AB + H A BPAcq BC + 2bq AB q AB = (25) 

The equations (16), (24) and (25) constitute the constraint equations. However, the equation (24) just confirms the 
equation (15) and then is not a independent constraint. The independent constraints are only (16) and (25). 
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IV. THE SYMMETRIC HYPERBOLICITY OF THE EVOLUTION EQUATIONS 



The principal parts of the evolution equations (17), (19) and (23) are, respectively 

F A = ^NC B C b c d {BC) ^ A + N CB d (BC) F A 

~Nc D r c d (BC) q DA + ^N^d {BC) r A 

_ foaADC a ~nN , .j.aADC a ~MN 
— * nN<J(DC)P + * MNO(dC) < 1 > 

and 

2 a/2 AT a/2 2N de 1 

= — ^ C 5(AC)PaS - -^Nif A D C^ C d {DC) p a B + g C, a Ed( A C)P&B + ^N CD ip A a d {DC) p aB 

~~2~4 (AC) ~ ~2i NvcA ^ d (ED)TB H g — d (AD) q CB 

Cl C DnN a ~ , tt CDMN a ~ 
= "AB 0(CD)PnN + LlAB 0(CD)QMN, 



where 



V2" 

* niV = — -iVC Cn £ Af + iV e n e AT ! 

x/2 

* mat — — ^-^s m? +—-l\tp €m cat , 

o o 

CDnAf 2 a/2 TV ng c w \/2 D>„C £ at , 2 ^ Fg D c w 1 at cd,. 71 N 

Was = g — £.4 es - "^"AtpA C £ s H g — </?f £a C b £ s + g-ZV (p^ £_b , 

and 

tt CDMN _ DM, AT , 

%/2 

Has = ~~24~ ~2T A 3 — ' 

The principal parts of the evolution equations can be expressed as a "matrix form" 

/ CaA \ ( \ / ( bE 

dt\ p aA = ^ aACD nN ^ACDMN q ~n< 
\q A B J V ®AB CD nN U AB CDMN J \ q M 

and then they are symmetric hyperbolic if the conditions 

QaACDnN _ -^nNCDaA 
qABCDuN _ -qjABCDnN 
U ABCDMN = y^MNCDAB 

are satisfied. Using the formulas 

JiA mi/- N /- A -NA -mbT NT A 

e - e Cm C& , e — e Cm U 

and 

— nAfCDaA _ a/2— n.D-^ aC _ Ajv + _CT>_ „_ AJV 

8 



we can find easily that the condition 



^aACDnN _ -^iNCDaA 



leads to 



Supposing 



and using 



— NQ aD (n c e NA _|_ ]\J CD e na £ NA _ V^J^ nD -^aC-^AN ^J^CD_ an _ AN 



N° D = N CD , 



rnb — mb 

t t , 



we find from (34) that 

e NA = e NA , N = ~N. 

Since 



and 



—ABC DM N = _V^ Ne MA^DC e BN + ^ N ^DM ^AC £ SiV ^ 



n.ABCDMN 2\/2N MD AC r BN V^2 AT AD MC r BN 2N EM D AC r BN , ^ \jCD AM 

" = 3 — V e e - -^j-^ e 6 3 — ¥>£ e e 3 ^ 



the condition 



leads to 



qABCDMN _ ^ABCDMN 



< W^N_^MD e AC e BN _^L Nip AD e MC e BN _ ™ EM ^ D ^ AC ^ BN + ^jyCD ^AM e BN 

3 24 3 3 



which means 



= -)fl N ^A^DC t BN + yl N ^DM e AC £ BN^ 



V AC = 3^ AC , 



and 

2N A B +N CD V A c^ BD = 0. 

Using (37) and (20) we obtain 



^ca^ AB = \ec B - 



And, (38) leads to trivial result = 0. 
Similarly, 



J^MNCDAB _ jy c MC c DA c NB _|_ jy ^MA^uCD Q N B + 2N DA ^ MC ^ NB 



24 24 

and 



jjABCDMAt = y^MNCDAB 



lead to e AM = -e MA , and N DC = N DC ', which are trivial. 

In summary, the conditions for the evolution equations of CaA, p aA and cjab together with the assumption (35) 
reduce to the reality conditions (36) on the metric and a condition (39) on the constant spatial vector i\) AB . 

In this case the polynomial Q p (p, q) in the equation (19) has the form 

Q P (p,q) 

64a ^ P PBC Ma P i B Ua P BP 

i N „i,aB~ C~ A 3JV , a B~ A~ C %N B Ac N aB ~ ~ AC 

+-^V PC QB PB QC -g^V V 1BC + ^ PBcq 

-^ aA p BC qBC - ^ AB PBcq aC + ^BCpF A q- B + ^ A WB C p BD ~q CD 
32cr 64cr 16(7 8cr 

/-aA~ ~BC N ~aA~ B N ~A ~aB 

~ 192^ C QBCq Wa q qB 48^ m 

~^r B ^C D r A q BD ^-r B ^ CA qc D q BD 
loa 64ct 

+d (B c)N CB p aA ^d {BC) NC, aB p^ A ^d (BC) Nr C q aA - ^d iBC) Nr B r A , (40) 
and the evolution equation (23) becomes 

Qab 

= -2V2NC aD ^D C d (AC) p aB - ^C D ipA C d (CD) p aB + 2N DE ^ E c CDd (AC) PaB + N CE VAD C aD d (CE) p aB 
-^-d {AC) ? B ^L M ™ d{DE) lf B + 2 -N CD d (AC) q DB 

+Q q (p,q) (41) 

where 



Qqip, q) 

25N , „ Dr „ 17N , 65N , ~ c ~ n , r~ n~ 

~~^T^ABP PDC ~ -TT-^ACjy PBD + -^-IpACPB PD + WD PA PCB 

b4a 64cr 64(7 

V2N AE c~ ~ed V2N de , c _ _ V2N DE , c _ _ y/2N DE c _ _ 
H 5 ipD PcbP H -. WA PcdPeb H -; Vt> PabPce H 5 PcbPae 

8(7 4(7 4(7 8(7 

AT „ DF „ AT „ „ r 35A^ r ^ AT „ r „ 2>hN „ „ r AT , r ^ FD ^ 

+ 7^ _e ^ B P <7l)£ + ^r^ B< ?C - -7^—Pc qAB + -TZ-PA qBC - -7^-PBCqA + -T-lpABlpE P qcD 

96(7 64(7 32(7 192(7 96(7 8(7 

+ -^->pAC^B D pD E q CE ~ -^ r i>Aci>DEPB E q <:D + ^ C ° VA E pD E q C B ~ — ° V A E PD B q C E 

QAcr 16cr 4cr cr 

V2A^ , „ q V2N DE , , _ c _ F v^A^ , c , f~ ~ 
H 3 WacWfePb qD ~ VacVfePe i «d ~ Vd Wa PEBqcF 

8(7 2(7 2(7 

V2A^ £ , c , f~ ~ V2A^ , c , „ p~ V2A^ , c , _ _ F 

H 5 WD WA PFBqCE H 5 VF VAEPD qCB 5 WF WAEPCBqD 

0(7 0(7 0(7 

HAT , „ ^ D A^ „ c _ D A^ „ AT , „ ^r, A^ „ c _ r> 

~TKn VABqCDH + T^-VACqB qD "I" tyACqBDy Trio V BCqADn + T^~VCDqB ■ qA 

192(7 16(7 48(7 192cr 48(7 

V^A^- , F ^ „ ViV2N CD , F ^ „ ^2^°-° , F ^ „ V2A^ AjD , F „ „ nr 

H — c ^4 qDBqcE vc qABqDE h 7, vc qoBqAE tta — <iEBq 

6(7 24(7 6(7 24(7 

+ || 9 (AC)^ B - 2V2d (AC) N V DC p DB - ^Id {DE) N VA D p E B + ^-d^NyAcv ^ B 

+ \d ( A C) N CD q DB - 2d [AC) N DE VD c p EB - d {DE) N ED VAcP <J B . (42) 
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V. CONCLUSIONS 



Now we have proved that the evolution equation (17), (19) and (41) constitute a symmetric hyperbolic system 
under the conditions (35), (36) and (39), and the constraint equations are (16) and (25). Here we suppose only the 
reality of the shift vector N CD and then obtain naturally the reality of the spatial metric as one of the conditions of 
the symmetric hyperbolic system rather than a independent one as in [23]. 

This work was supported by the National Science Foundation of China No. 10175032. 
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